Non-equilibrium dynamics in an interacting nanoparticle system 
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Non-equilibrium dynamics in an interacting Fe-C nanoparticle sample, exhibiting a low tem- 
perature spin glass like phase, has been studied by low frequency ac-susceptibility and magnetic 
relaxation experiments. The non-equilibrium behavior shows characteristic spin glass features, but 
some qualitative differences exist. The nature of these differences is discussed. 

I. INTRODUCTION 

The dynamics of interacting nanoparticle systems has been subject of considerable interest concerning the existence 
of a low temperature spin glass phase. Evidence for a phase transition from a high temperature superparamagnetic 
phase to a low temperature spin glass like_phase has been given from reports on critical slowing dowru and a divergent 
behavior of the nonlinear susceptibility.Ha 

Spin glasses have been intensively investigated during the two last decades, not only concerning the phase transition 
but also the nature of the spin glass phase. The non-equilibrium dynamics in the spin glass phase has been exten- 
sively studied by conventional and also more sophisticated dc-magnetic relaxation and low frequency ac-susceptibility 
experiments.l3 A non-equilibrium character of the low temperature spi]i,|glass like phase of interacting particle sys- 
tems has beefl revealed from measurements of dc-magnetic relaxatioroB and relaxation of the low frequency ac- 
susceptibility.H In this paper, the non-equilibrium dynamics in the spin glass like phase is furtherj-jekLcidated and 
special interest is put into a dynamic memory effect, which has recently been observed in spin glasses.Ll~t2l In a "mem- 
ory experiment" the ac-susceptibility is measured at a low frequency, and the cooling of the spin glass sample is halted 
at one (or more) temperatures. During a halt, x slowly decays, but when cooling is resumed x gradually regains the 
amplitude of a continuous cooling experiment. Upon the following heating, x shows dips at the temperatures where 
the temporary halts occurred; a memory of the cooling history has been imprinted in the spin structure. (For an 
illustration see Fig. 11a Ref. ^.) A similar memory of the cooling process has recently been reported in an interacting 
nano-particle system.Eil 

II. THEORETICAL BACKGROUND 



For spin glasses, the nonequilibrium dynamics has heen interpreted within a phenomenological real space model,El 
adopting important concepts from the droplet model.ll3 We will use the same real space model to discuss similarities 
and dissimilarities between the low temperature phase of interacting particle systems and that of spin glasses. 

The droplet model was derived for a short-range Ising spin glass, but important concepts like domain growth, chaos 
with temperature, and overlap length should be applicable also for particle systems exhibiting strong dipole-dipole 
interaction and random orientations of the anisotropy axes. Chaos with temperature means that a small temperature 
shift changes the equilibrium configuration of the magnetic moments completely on sufficiently long length scales. The 
length scale, up to which no essential change in configuration of the equilibrium state is observed after a temperature 
step AT, is called the overlap length 1{AT). The development towards equilibrium is governed by the growth of 
equilibrium domains. The typical domain size after a time iw at a constant temperature T is 



A(r) 

where r is the relaxation time of an individual magnetic moment, A(T) sets the free energy scale and "0 is a barrier 
exponent. For spin glasses, the atomic relaxation time is of the order of 10~^^ s independent of temperature, while 
for magnetic particles, the individual particle relaxation time is given by an Arrhenius law as 
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where tq ~ lO"""^^ — 10~^ s, K is an anisotropy constant, V the particle volume and thus KV the anisotropy energy 
barrier. Due to the inevitable polydispersivity of particle systems the individual particle relaxation time will also be 
distributed, and due to the exponential factor this distribution will be broader than the distribution of anisotropy 
energy barriers. 

In the droplet model there exists only two degenerate equilibrium spin configurations, ^! and its spin reversal 
counterpart ^ . It is thus possible to map all spins to either of the two equilibrium configurations. Let us now consider 
a spin glass quenched to a temperature Ti lower than the spin glass transition temperature Tg. At t = the spins will 
randomly belong to either of the two equilibrium states, resulting in fractal domains of many sizes. The subsequent 
equilibration process at t > is governed by droplet excitationsyielding domain growth to a typical length scale, which 
depends on temperature and wait time <wi according to Eq.(|]). After this time, fractal structures typically smaller 
than R{Ti, t^^) have become equilibrated while structures on longer length scales persist. If the system thereafter is 
quenched to a lower temperature T2, the spins can be mapped to the equilibrium configuration at T2 yielding a new 
fractal domain structure. The domain structure achieved at temperature Ti still fits to the equilibrium configuration 
at T2 on length scales smaller than the overlap length l{Ti — T2). The domain growth at T2 will start at the overlap 
length, increase with wait time iwj: a-nd end up at the size i?(r2,tw2)- Heating the system back to Ti, the fractal 
domain growth that occurred at T2 has only introduced a new and dispersed domain structure at Ti on length scales 
l{Ti — T2) < R < i?(T2,tw2). Note that the large length scale structures, R > R{Ti,tv,^), essentially persist. 

The time-dependent response m(iobs) to a weak magnetic field applied at tobs = 0, is due to a continuous magneti- 
zation process, governed by polarization of droplets of size 

i(T.,...,.fIlMV'*. (3) 



A(T) 

Since L(T, iobs) grows with the same logarithmic rate as i?(T, t^), the relevant droplet excitations and the actual 
domain sizes become comparably large at time scales In tobs ~ In tw . For In tobs ^ In tw the relevant excitations occur 
mainly within equilibrated regions, while for Iniobs 3> Ini^ these excitations occur on length scales of the order of the 
growing domain size, and involves domain walls, yielding a nonequilibrium response. A crossover from equilibrium 
to nonequilibrium dynamics occurs for In tobs ~ Intw, seen as a maximum (or only a bump) in the relaxation rate 
S{t) = h~^dm{t)/dlnt vs. Int curves. 

Dc-relaxatiofl, and ac-susceptibility experiments are related through the relations: M{t) w x'(^) ^^'^ S{t) cx 
when t = l/tj.E£l Hence, the aging behavior is also observed in low frequency ac-susceptibility measurements. Different 
from dc-relaxation measurements is that the observation time is constant, tobs = 1/^: implying that the probing length 
scale L{T, 1/lo) is fixed for a given temperature. 



III. EXPERIMENTAL 



The sample consisted of ferromagnetic particles of the amorphous alloy Fei-xCx {x «0.22) prepared by thermal 
decomposition of Fe(C0)5 in an organic liquid (decahn) in the presence of surfactant molecules (oleic acid) as described 



in Ref. 14. After the preparation, the carrier liquid was evaporated and the particles were transferred to xylene in 
an oxygen-free environment resulting in a ferrofluid with a particle concentration of 5 vol% (determined by atomic 
absorption spectroscopy) . The particles were separated by the surfactant coating and could only interact via magnetic 
dipole-dipole interactions. During measurements, the ferrofluid was contained in a small sapphire cup sealed with 
epoxy glue to prevent oxidation of the particles. The sample was only measured and exposed to magnetic fields well 
below the melting point of xylene (« 248 K). A droplet of the ferrofluid was dripped onto a grid for transmission 
electron microscopy and was left to oxidize prior to the study. The electron micrographs revealed a nearly spherical 
particle shape. The volume-weighted size distribution (after correction for the change in density due to oxidation) was 
described well by the log-normal distribution, f{V)dV = {V2T:avV)~^ exp[— In^ (V/Vm) / {'2(Ty)]dV , with the median 
volume Vm — 8.6 • 10^^^ m'^ and logarithmic standard deviation cry — 0.19. These parameters correspond to a 
log-normal volume- weighted distribution of particle diameters with the median value d„i — {GYm/nY^^ ~ 5.5 nm and 
ad — cry/3 =_(1.QE2. This particle size is slightly larger and the value of cry is slightly lower than those reported in 
recent studiesElaa performed on a different batch of particles. The 5 vol% sample has been subject to a detailed study 
of the correlations and magnetic dynamics using a wide range of observation times and temperatures Collective 
behavior is observed at temperatures below 40 K. j— , 

A non-commercial low-field SQUID magnetometeiO was used for the measurements. All ac-susceptibility measure- 
ments were performed at a frequency of 510 mHz and a rms value of the ac-field of 0.1 Oe. The cooling and heating 
rate was 0.25 K/min. For the relaxation measurements the dc field was 0.05 Oe, while the background field was less 
than 1 mOe. 
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IV. RESULTS AND DISCUSSION 



A. General behavior 

Fig. 1^ shows the ac susceptibihty of the particle sample measured while cooling. The onset of the out-of-phase 
component of the susceptibility is less sharp than for spin glasses. Nevertheless, there exists collective spin glass 
like dynamics in the particle system, which can be studied by magnetic aging experiments. The 'conventional' aging 
experiment is a measurement of the wait time dependence of the response to a field change, after either field cooling 
(thermorcmanent magnetization, TRM) or zero field cooling (ZFC). In Fig. ^ the relaxation rate from two ZFC 
magnetic relaxation measurements is shown. The particle system was directly cooled from a high temperature to 
30 K where it was kept a wait time of 300 and 3000 s, respectively, before the magnetic field was applied. A clear 
wait time dependence is seen; the relaxation rate displays a maximum at Iniobs ~ lnt„. Both particle systems and 
spin glasses exhibit magnetic aging that affects the response function in a similar way, but the wait time dependence 
appears weaker for the particle system. Comparing the ratio between the magnitude of S{t) at the maximum and 
at the short time limit (0.3 s in our current experiment) wc find a value of order 1.3 in the current sample, whereas 
most spin glasses would show a value > 2, at a corresponding temperature and wait time. The illustrated differences 
between the dynamics of a particle system and an archetypical spin glass can be assigned to the wide distribution of 
particle relaxation times and that some particles may relax independently of the collective spin glass phase due to 
sample inhomogeneities. 

The relaxation in the ac-susceptibility, after cooling the sample to 23 K with the same cooling rate as in the reference 
measurement, is shown in Fig. ^. The relaxation in absolute units is larger for x' than for x" and, in contrast to spin 
glasseall where the aging phenomena are best exposed in x", the effects of aging are more clearly seen in x' for this 
particle system. Measuring the ac-susceptibility with slower cooling rates gives lower values of the ac-susceptibility. 
The effect of a slower cooling rate is largest just below the transition temperature where the effect of aging also is 
largest. 

B. Temperature cycling 

One expedmental procedure that has been used for spin glasses to confirm the overlap length concept is tempera- 
ture cyclingJij The measurements are similar to conventional aging experiments, but after the wait time i^i at the 
measuring temperature Tm, a temperature change AT is made, and the system is exposed to a second wait time 
before changing the temperature back to Tm, where the field is applied and the magnetization is recorded as a 
function of time. 

Fig. ^ shows results from measurements with negative temperature cycling, at Tm = 30 K, using twi= 3000 and 
=10000 s. For > AT > —2 K the additional wait time at the cycling temperature mainly shifts the maximum in 
the relaxation rate to longer times. The shift is largest for the smallest AT and the maximum returns continuously 
towards iwi as AT is increased. For AT = — 2 K, the maximum appears at Ini « Iniwi and the relaxation rate behaves 
rather similar to the curve without the cycling only showing a somewhat larger magnitude (see Fig For a larger 
temperature step, AT = —4 K, the relaxation rate is enhanced at times shorter than , while the relaxation appears 
unaffected at longer time scales. The behavior can be interpreted in terms of an interplay between the domain growth 
and the overlap length: for small temperature steps, the domain size i?i(Ti„,twi) attained during t^-^ at T^ is shorter 
than the overlap length at the cycling temperature and the domain growth proceeds essentially unaffected by the 
temperature change, only with a slower rate. For sufficiently large temperature steps the overlap length is shorter than 
i?i(Tm,iwi) and the domain growth at the lower temperature creates a new domain structure on short length scales, 
^2(Tm + AT, twg)- For the temperatures and wait times used in the experiments i?2(T„i + AT, ^^2) < -Ri(?m,iwi)- 
Returning to Tm, the new R2 structures do not overlap with the equilibrium configuration and yield the apparent 
non-equilibrium nature of the dynamics on short time scales (cf. Fig. ^ for AT = —4 K). 

Measurements with positive temperature cycling with Tm = 27 K are shown in Fig. |^. For small temperature steps 
(AT ^ 2 K) the maximum of the relaxation rate remains at In i « In t^i , but for larger temperature steps it is seen 
that the relaxation rate approaches the relaxation rate of an aging measurement with short wait time. This can again 
be explained by the overlap length becoming shorter than i?,i(Tm, iwi) when AT ^ 2 K. Since the domain growth rate 
increases fast with increasing temperature, the system will then look more and more reinitialized when returning to 
Tm after larger AT and/or longer tw2- 
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C. Memory 



The experimental procedure for a memory experiment is illustrated in Fig. ^. x'C^) is recorded during cooling, 
employing different halts at constant temperature, followed by continuous heating. Fig. Q shows results from three 
different memory experiments: the first with a single temporary halt at 33 K for 1 h and 30 minutes, the second 
with a single temporary halt at 23 K for 10 h, and the third with temporary halts at both 33 K and 23 K. The 
exposed curves show the difference between the measured and a reference curve. Reference curves are taken from a 
measurement with continuous cooling followed by continuous heating. Similar to spin glasses, a memory of the cooling 
history is observed on heating for both one and two temporary halts. It is seen that on the low temperature side, the 
susceptibility approaches the reference level slower than on the high temperature side. This asymmetry of the dip is 
more pronounced in the particle system than for spin glasses. The relaxation at 23 K after a temporary halt at 33 K 
is equally large as the relaxation at 23 K without earlier halt, except that it starts at a lower level. In fact, x' ~ Xrof > 
measured on heating for the experiment with two temporary halts at 33 K and 23 K, is just the sum of x' ~ xJcf '^'^ 
the two experiments with a single temporary halt at 33 K and 23 K, respectively. 

In Fig. 1^ is shown a memory experiment similar to the one in Fig. ^ except that the two temporary halts are closer 
in temperature. The ac-susceptibility has been measured with temporary halts on cooling at 33 K for lh30min and at 
28 K for 7 h. On heating, the double halt experiment only shows one dip. However, this dip (x' — Xrof ) still constitutes 
the sum of the two heating curves with one halt, as was the case for the experiments shown in Fig. ^. Also for this 
experiment, the ac-relaxation curve at the lower temperature is the same with and without temporary halt at 33 K, 
except for a constant. 

We have performed supplementary experiments to the ones shown in Fig. |^ and ^ by measuring the dc-relaxation. 
In Fig. ^, the relaxation rate of three dc-relaxation curves measured at 33 K are shown. In the three experiments 
the system was quenched to 33 K and then aged for IhSOmin. In one of the measurements the relaxation was 
measured immediately after this wait time, but in the two other measurements a negative temperature cycling was 
performed to 23 K for 10 h and to 28 K for 7 h, respectively. The dynamics is sizeably affected by the temperature 
cycling to 28 K, but only weakly affected by the temperature cycling to 23 K. The observed increase of the relaxation 
rate indicates that the equili brium structure created at 33 K is erased on short length scales after the temperature 



cycling, as discussed in Sec. IV B. On time scales shorter than 10 s, S'(i) is even lower for the experiment with a 
temperature cycling to 28 K than for the other two experiments. We conclude that there is an overlap between the 
equilibrium configuration of the magnetic moments at 33 K and 28 K on length scales shorter than L(33 K, 10 s). 
Let us now compare these dc-measurements with the corresponding ac-measurements. The in-phase component of 
the ac-susceptibility at a frequency of 510 mHz records the integrated response corresponding to observation times 
1/w « 0.3 s and shorter. In the ac-experiment with two temporary halts at 33 and at 23 K (Fig. |^), no effect of the 
aging at 23 K is seen when heating through 33 K. However, x'(33 K) measured on heating for the ac- measurement 
with temporary halts at 33 and 28 K is lower than x'(33 K) for the measurement with only a single halt at 33 K. A 
lower level of x'(^) indicates a more equilibrated system and thereby also a lower relaxation rate. This supports the 
above conclusion of some overlap between the equilibrium configurations at 28 K and 33 K on the length scales seen 
by observation times shorter than 0.3 s. 

Fig. |l^ shows the relaxation rate measured at 28 and 23 K after a quench to Tm followed by a wait time of 7 h and 
10 h, respectively. The same measurements except for a temporary halt at 33 K for lh30min are also shown. The 
relaxation rate is lower for the measurement with a temporary halt at 33 K for both temperatures, even though the 
difference between a temporary halt and a direct quench is larger at 28 K. These measurements, as the corresponding 
ac-measurements, show that there is some overlap on short length scales between the equilibrium configuration at 
33 K and all lower temperatures for which we have measured the ac-susceptibility. 



V. CONCLUSIONS 



The non-equilibrium dynamics in the low temperature spin glass like phase of an interacting Fe-C nano-particle 
sample is found to largely mimic the corresponding spin glass dynamics. The observed differences may be accounted 
for by the strongly temperature dependent and widely distributed relaxation times of the particle magnetic moments 
compared to the temperature independent and monodispersed relaxation times of the spins in a spin glass. Within a 
droplet scaling picture of a particle system, these factors strongly affects the associated length scales and growth rates 
of the domains and droplet excitations as well as the roughness of the domain walls. In addition, there might exist 
individual particles that do not take part in the collective low temperature spin glass phase but relax independently. 
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FIG. 1. The ac-susceptibility vs. temperature at / = 510 mHz. 
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FIG. 2. The relaxation rate of the zero field cooled magnetization vs. time at T = 30 K. The wait time is 300 s (circles) and 
3000 s (squares). 
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FIG. 3. Ac-susceptibility vs. time at T = 23 K. Same units as in Fig 0, / = 510 mHz. 
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FIG. 4. Relaxation rate S{€) vs. time. The sample has been aged twi = 3000 s at Tm = 30 K and afterward subjected to a 
negative temperature cycling AT during = 10000 s, immediately prior to the application of the field H — 0.05 Oe; Circles: 
AT = 0; Pluses AT = -2 K; Triangles: AT = -4 K. 
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FIG. 5. Relaxation rate S{t) vs. time. The sample has been aged twi = 3000 s at Tm = 27 K and afterward subjected to 
a positive temperature cycling AT during a time t^^ i immediately prior to the application of the field H = 0.05 Oe; Circles: 
Ar = 0; Pluses AT = 2 K, = 5 s; Triangles: AT = 4 K, = 5 s; Diamonds: AT = 6 K, twa = 30 s. Squares are data 
obtained after the sample has been aged at 27 K for 5 s only. 
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FIG. 6. The experimental procedure in an ac-susceptibility memory experiment with two temporary halts during cooling 




FIG. 7. x'{T) — x'ici{T) vs. T measured while cooling (solid symbols) or heating (open symbols). Circles, the cooling was 
halted at 33 K for 1 h 30 min; squares, the cooling was halted at 23 K for lOh and triangles, the cooliiig was halted at Ti = 33 K 
for twi = 1 h 30 min and at 72 = 23 K for = 10 h. Units and frequency are the same as in Fig. hi 
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FIG. 8. x'{T) — Xrof(^) vs. T measured while cooling (solid symbols) or heating (open symbols). Circles, the cooling was 
halted at 33 K for lh30min; squares, the cooling was halted at 28 K for 7 h and triangles, the cooling was halted at Ti = 33 K 
for twi = 1 h 30 min and at T2 = 28 K for twi = 7 h. Units and frequency are the same as in Fig. ^ 



10.5 




FIG. 9. Relaxation rate S{t) vs. log(t). The sample has been aged for lh30min at Tm = 33 K and afterward subjected to a 
negative temperature cycling AT during t^^ . Circles: AT = 0; Pluses: T^ + AT = 28 K, t^^ = 7 h; Triangles: T^^ + AT = 23 K, 
tw, = 10 h. H = 0.05 Oe. 
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FIG. 10. Relaxation rate S{t) vs. log(t). The sample has been aged, after a quench, for = 7 h at Tm = 28 K (circles) and 
after a quench with a temporary halt at 33 K for lh30min (triangles), H = 0.05 Oe. Inset: Tm = 23 K and tw ~ 10 h. 
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